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1. 


INTRODUCTION 


This  report  evaluates  the  surface  evaporative  duct  height 
field  produced  by  the  global  surface  contact  layer  interface 
model  at  the  Fleet  Numerical  Oceanography  Center  (FNOC)  .  This 
evaluation  is  intended  to  inform  users  of  the  product's  ability 
to  provide  a  representative  analyses  and  forecast (s)  of  global 
surface  ducting  conditions . 

Analysis  and  forecast  fields  from  the  Navy  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS)  are  required  as 
input  for  the  creation  of  the  surface  evaporative  duct  height 
field.  This  evaluation  was  conducted  prior  to  the  implementation 
of  NOGAPS  3.0  by  the  FLENUMOCEANCEN  in  January  1988.  Thus,  some 
of  the  conclusions  discussed  below  may  no  longer  be  valid.  An 
analysis  similar  to  the  one  described  below  could  be  performed  in 
order  to  assess  the  impact  of  the  new  global  prediction  system. 

A  surface  evaporative  duct  is  said  to  exist  when  a  rapid 
decrease  in  specific  humidity  q  near  the  surface  results  in  a 
decrease  in  the  modified  index  of  refraction  M  with  altitude 
(Richter  and  Hitney,  1979) .  The  relationship  between  M  and  q  is 

N  =  77.6  P/T  +  3.73  X  106  e/T2 

M  =  N  +  z/a 
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where 

q  =  0 . 622e/  (P  -  0.378e) 

and  where  P  is  the  pressure  (mb) ,  T  is  the  absolute  temperature 
(K)  ,  e  is  the  vapor  pressure  (mb) ,  z  is  the  altitude  (m) ,  and  a 
is  the  radius  of  the  earth. 
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When  dM/dz  <  0,  microwave  energy  may  be  trapped  near  the 
surface  and  radar  return  strength  may  be  enhanced.  Microwave 


rays  within  the  duct,  traveling  in  directions  near  grazing 
incidence  with  the  duct  boundaries,  may  be  unable  to  penetrate 
the  duct  and  may  be  reflected  by  the  duct  boundary  and  the  sea 
surface.  Thus,  such  rays  will  be  able  to  travel  long  distances 
downrange  with  minimal  loss  in  signal  strength.  The  critical 
incidence  angle  for  trapping  depends  on  frequency,  duct  height, 
and  duct  strength  (Ko  et  al . ,  1983).  Low  frequencies,  in 
general,  do  not  have  as  large  a  variation  in  extended  range 
with  varying  duct  heights  as  higher  frequencies,  and  higher 
frequencies  do  not  require  as  large  a  duct  thickness  for  the 
same  range  enhancement  as  lower  frequencies  (Polios  and  Hitney, 
1979) .  The  height  of  the  duct  is  where  the  lowest  inflection 
appears  in  the  M  profile  and  the  strength  of  the  duct  is  the 
magnitude  of  dM/dz  (Ko  et  al . ,  1983) .  Sea  clutter  is  also 
enhanced  in  ducting  situations. 

Many  circumstances  which  would  lead  to  large  gradients  in 
q  near  the  surface  of  the  ocean  can  be  envisioned.  In  fact, 
surface  evaporative  ducting  is  nearly  always  present  over  almost 
all  oceans  (Richter  and  Hitney,  1979) . 

The  existence  of  a  surface  evaporative  duct  has  many  serious 
operational  implications.  One  of  the  most  serious  and  obvious 
implications  is  that  the  extended  range  of  the  radar  will  allow 
distant  hostile  units  to  detect  the  radar  (and  to  be  detected  by 
the  radar) .  Thus,  the  existence  of  a  surface  evaporative  duct 
is  an  important  component  of  the  decision  on  whether  or  not  to 
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radiate.  Another  consequence  of  the  duct  is  that  a  radar  "hole" 
will  exist  between  the  top  of  the  duct  and  the  base  of  the  lobes 
unaffected  by  the  duct  (see  Figure  1) .  Hostile  aircraft  or 
missiles  will  thus  be  able  to  approach  the  radar  undetected. 
Conversely,  knowledge  of  the  refractivity  profile  in  the 
vicinity  of  enemy  radars  can  be  turned  to  tactical  advantage. 


Figure  1.  Schematic  of  shipboard  radar  coverage  under  normal 
atmospheric  conditions,  no  ducting.  A;  and  changes  in  the 
presence  of  an  evaporation  duct,  B.  under  ducting  conditions 
an  aircraft  in  the  radar  "hole"  can  not  be  seen  by  the  ship's 
radar,  as  indicated  in  B. 


Enhanced  sea  clutter  may  also  pose  a  threat .  Sea  clutter 
may  have  a  cross-section  larger  than  a  low-flying  aircraft  or 
missile,  thereby  hiding  it  from  the  radar.  Also,  automated  tracking 
systems  can  be  saturated  and  disabled  by  elements  of  sea  clutter. 
(Snyder,  1979) .  An  unusual  amount  of  sea  clutter,  however,  will 
provide  the  operator  with  evidence  of  ducting  conditions.  This 
may  be  the  only  indication  of  ducting  if  the  refractivity  profile 
is  not  known. 

There  are  several  sources  of  operational  information  on 
refractivity.  The  easiest  to  use  is  the  refractive  effects 
guidebook  (REG)  in  which  synoptic  typing  is  used  to  associate 
certain  refractivity  profiles  with  climatological  patterns. 

Surface  evaporative  ducts,  however,  are  beyond  the  scope  of 
the  Refractive  Effects  Guidebook  (REG)  published  by  NEPRF 
(Rosenthal,  1976) .  Another  source  of  information  is  the 
Integrated  Refractive  Effects  Prediction  System  (IREPS) ,  which 
requires  as  input  a  radiosonde  sounding.  When  an  elevated 
trapping  layer  overlies  a  surface  evaporative  duct,  the  behavior 
of  the  radar  will  be  dramatically  effected  (Sketchly,  1979) .  In 
such  a  situation,  however,  IREPS  will  assume  that  the  surface 
duct  is  the  dominant  feature  (Paulus  and  Hitney,  1979) .  Finally, 
the  sea  clutter  algorithms  in  IREPS  are  more  consistent  with 
ducts  based  at  the  surface  formed  by  elevated  refractive  layers 
(i.e.,  a  duct  in  which  the  elevated  minimum  in  the  M  profile  is 
smaller  than  the  surface  value]  (Snyder,  1979) .  Although  IREPS 
is  capable  of  diagnosing  the  thickness  of  the  surface  evaporative 
duct,  the  only  product  specifically  designed  for  the  prediction 
of  surface  evaporative  ducts  is  the  surface  evaporative  duct 
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height  field  produced  at  the  FNOC  by  the  Global  Surface  Contact 
Layer  Interface  (GSCLI)  model. 


2.  GLOBAL  SURFACE  CONTACT  LAYER  INTERFACE  MODEL 

The  GSCLI  model  is  a  generalized  similarity  theory  model 
for  the  entire  planetary  boundary  layer  (PBL)  in  which  three 
nondimensional  parameters  are  assumed  to  completely  determine 
the  momentum  and  mass  fields  within  the  PBL.  The  first  two 
parameters  are  the  surface  Rossby  number  Ro  and  the  temperature 
stratification  parameter  S: 

Ro  =  G(fz0)_1 
S  =  gA©v  (©vfG)_1 

where 

Aev  =  ev  -  ev 

X  o 

Aq  =  qx  -  qQ 


and  where  G  is  the  500  m  gradient  wind,  f  is  the  Coriolis  param¬ 
eter,  zQ  is  the  roughness  length  (which  is  a  function  of  wind 
speed),  and  g  is  the  acceleration  due  to  gravity.  The  subscript 
”x"  refers  to  the  top  of  the  surface  layer  at  h=BG/f  (B  is  a 
constant)  while  the  subscript  "o"  refers  to  the  surface.  Ro 
and  S  determine  the  barotropic  PBL.  The  third  parameter  V 
introduces  baroclinicity : 


where 

t 

. 

'  and 
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M  =  H'/L 
H  '  =  ku*(f)_1 

ev  u*3 

L  =  - 

kg  (u*0*  +  0.609  u*q*) 


5 


and  where  k  is  von  Karmen' s  constant.  The  scale  velocity  u*, 
scale  virtual  potential  temperature  ©*,  and  scale  moisture  q* 
are  qiven  by 

u*  =  CdG 

0*  =  u*2  9v(k2gL)_1 

q*  =  Cp©v  (  ACp  ©v/AAq) 

where  CD  is  the  drag  coefficient  (which  is  a  function  of  z/L) , 

C_  is  the  specific  heat  of  air  at  constant  pressure,  and  A  is  the 

r 

latent  heat  of  vaporization.  L  is  the  Monin-Obukhov  length  scale 
which  is  a  ratio  of  mechanical  to  bouyant  turbulent  kinetic 
energy  production.  Ro,  S,  and  m  are  determined  from  output  from 
the  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS) .  Note  that  the  GSCLI  is  purely  a  diagnostic  system; 
forecast  fields  from  the  GSCLI  are  generated  by  using  NOGAPS 
forecast  fields  for  determination  of  Ro,  S,  and  y.  Within  the 
surface  layer,  the  scaling  quantities  are  related  to  the  mean 
gradients  by 


where  the  diabatic  flux-profile  function  4>  is  given  by 

z 

<J)  =  0.74  +  4.7  (-)  stable 

L 

z  -1/2 

=  0.74  (1-9  (-)  )  unstable 

L 

and  where  z  is  a  reference  height .  Given  the  scaling  quantities 
q*  and  0*,  an  expression  for  the  surface  evaporative  duct  height 
(derived  in  Appendix  A)  can  be  obtained: 

C  q*  <£(z/L)  -  T*  <p(z/L) 

5  =  - ^ - 

b2  TZ/AP  +  0.243 

where  C  and  A  are  constants  and  P  is  the  pressure  (mb) .  b2  is 
the  critical  vertical  gradient  in  N  for  trapping  (=  -157  N/Km) . 

5  is  produced  operationally  over  all  ocean  regions  at  6  hour 
intervals  from  tau  0  to  tau  36.  The  record  name  is  A90 .  Several 
related  fields  are  also  produced,  including  A92,  the  surface  N 
value  minus  the  value  at  the  level  where  dN/dz  =  -157  N/Km. 

A  discussion  of  surface  layer  similarity  theory  as  it 
pertains  to  near-surface  gradients  in  N  appears  in  Appendix  B. 

The  primary  NOGAPS -derived  quantities  which  influence  the 
GSCLI  are  the  low-level  pressure  gradient  and,  in  the  tropics, 
the  PBL  wind  field  (which  determines  G) ,  the  sea-surface  temper¬ 
ature  field,  and  the  850  mb  temperature  field.  Thus,  the  amount 
of  "synoptic"  information  available  to  the  GSCLI  is  limited. 
Determination  of  atmospheric  stability  in  the  GSCLI  relies  on  an 
extrapolation  of  the  850  mb  temperature  to  1000  mb  (using  the 
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the  1000-850  mb  thickness)  and  a  further  extrapolation  of  the 
difference  between  the  extrapolated  1000  mb  temperature  and  the 
sea  surface  temperature  to  arrive  at  the  air/sea  temperature 
difference.  Stability  is  critical  in  determination  of  6 .  It 
plays  a  role  in  determining  L,  in  selecting  the  appropriate  <£, 
and  in  computing  q*  and  9*.  Several  constraints  on  L,  z/L,  and 
0  are  required  by  the  formulation  of  the  GSCLI  model  which  limit 
its  generality.  A  number  of  ad  hoc  corrections  have  also  been 
added  in  order  to  improve  model  performance  in  certain  synoptic 
situations  (e.g.,  cold  air  outbreaks). 

In  light  of  uncertainties  concerning  the  quality  of  the  model 
products,  a  study  has  been  undertaken  at  the  Naval  Environmental 
Prediction  Research  Facility  (NEPRF)  in  order  to  evaluate  the 
response  of  the  A90  product  to  changes  in  the  synoptic  pattern. 

Although  the  analysis  (tau=0)  field  is  subject  to  the  same 
uncertainties  as  the  forecast  fields,  the  lack  of  any  independent 
measurements  of  surface  evaporative  ducting  necessitates  verifi¬ 
cation  of  the  forecast  against  the  analyses . 

The  analyses  presented  in  Appendices  A  and  B  of  the  expres¬ 
sion  used  for  surface  evaporative  duct  height  in  the  GSCLI  model 
clearly  demonstrate  a  serious  shortcoming  of  this  expression. 


Specifically,  when  the  ambient  conditions  are  such  that  dN/dz 
never  approaches  the  critical  value,  the  GSCLI  model  falsely 
produces  reasonable  appearing  duct  heights.  When  dN/dz  is  much 
greater  than  the  critical  value,  the  duct  height  should  be  much 
larger  than  the  value  produced  by  the  GSCLI  model .  When  dN/dz 
is  much  less  than  the  critical  value,  no  duct  exists.  Based 
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on  the  development  in  Appendix  B,  it  is  strongly  recommended 
that  the  similarity  theory  approach  for  potential  refractivity 
be  adopted  under  stable  conditions  for  the  operational  surface 
evaporative  duct  height  product . 

3.  EVALUATION  OF  THE  SURFACE  EVAPORATIVE  DUCT  HEIGHT  PRODUCT 
During  a  period  of  several  months  over  three  seasons,  root 
mean  square  (RMS)  errors  were  computed  for  the  A90  field  each 
day  at  both  00Z  and  12Z  for  the  12,  24,  and  36  hour  forecasts. 
The  error  fields  were  constructed  by  subtraction  of  12,  24,  and 
36  hour  forecasts  valid  at  the  analysis  time  from  the  analysis 
field.  In  addition,  the  maximum  error,  location  of  the  maximum 
error,  and  higher  moments  of  the  error  distribution  (variance, 
skewness,  and  kurtosis)  were  computed.  The  RMS  errors  were 
relatively  insensitive  to  the  length  of  the  forecast;  RMS  errors 
in  the  12  hour  forecast  rarely  departed  very  far  from  the  range 
of  2. 9-3.1  m  with  an  increase  of  only  about  15%  by  36  hours. 
Error  variance  tended  to  be  in  the  range  from  9  to  11  m  with  an 
increase  of  2-3  m  by  36  hours.  Maximum  errors  rarely  departed 
very  far  from  the  range  of  17-23  m  with  no  discernible  trend  for 
increasing  length  of  forecast.  No  diurnal  or  seasonal  variation 
was  detectable. 

3.1  Characteristic  Features 

Surface  evaporative  ducts  (SED)  are  primarily  the  result 
of  large  near  surface  gradients  in  specific  humidity  (or  vapor 
pressure) .  Consequently,  one  might  expect  that  the  most  active 


surface  evaporative  ducting  regions  would  be  in  the  tropics 
where  high  temperatures  and  relative  humidities  lead  to  high 
specific  humidities.  While  SED  is  by  no  means  confined  to  the 
tropics,  the  greatest  frequency  of  occurrence  (in  the  forecast 
and  analyzed  A90  field)  is  in  the  region  from  30  degrees  north 
to  30  degrees  south. 

Tropical  cyclones,  in  particular,  produce  a  well-defined 
signature  in  the  SED  field.  Typically,  the  tropical  cyclone  is 
associated  with  a  local  maximum  in  duct  height.  High  wind  speeds 
promote  mechanical  mixing  which  results  in  a  well  mixed  PBL  with 
large  gradients  in  specific  humidity  and  virtual  potential 
temperature  confined  to  the  capping  inversion  at  the  PBL  top. 
Large  errors  in  SED  forecast  fields  are  associated  with  tropical 
cyclones.  The  tropical  cyclone  bogus  is  present  in  the  analysis 
but  not  present  in  the  forecast  fields .  This  leads  to  large 
differences  between  analysis  and  forecast  fields  in  the  vicinity 
of  a  tropical  cyclone. 

A  characteristic  feature  associated  with  the  importance  of 
specific  humidity  is  local  maxima  in  duct  height  due  to  ocean 
currents.  This  is  particularly  true  of  warm  currents  in  the 
winter  hemisphere.  The  relatively  warm  water  promotes  large 
surface  moisture  fluxes  so  that  the  air  near  the  surface  is  quite 
moist  while  the  air  above  is  cool  and  dry.  The  gulf  stream,  in 
particular,  is  quite  evident  in  the  SED  fields  in  January;  the 
maximum  SED  values  coincide  with  the  maximum  in  the  NOGAPS  sea 
surface  temperature  (SST)  field. 
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A  feature  conspicuous  by  its  absence  is  a  local  maximum  in 
SED  height  and  frequency  off  the  western  coasts  of  continents  in 
the  subtropical  latitudes  of  the  summer  hemisphere.  For  example, 
surface  ducts  are  quite  commonly  observed  in  the  summer  in  the 
southern  California  operating  area  (Paulus  and  Hitney,  1979) ,  but 
do  not  appear  in  the  analyzed  or  forecast  SED  fields. 

3.2  Midlatitudes 

Midlatitude  response  of  the  SED  product  to  storms  traversing 
oceanic  areas  exhibits  a  pronounced  geographic  variability.  In 
the  North  Atlantic  and  South  Pacific,  the  product  responds  well 
to  oceanic  extratropical  cyclones  while  it  is  deficient  in  the 
North  Pacific  and  South  Atlantic  oceans.  There  is  no  apparent 
diurnal  difference  in  forecast  quality.  Examples  of  this 
geographic  bias  are  presented  below.  Possible  causes  of  the 
geographical  bias  are  discussed  in  the  conclufion. 

Examination  of  a  number  of  cases  reveals  several  general 
features.  Over  the  North  Atlantic  ocean,  the  forecast  ducting 
response  is  typically  on  the  south  east  side  of  the  low  pressure 
area  (possibly  in  the  vicinity  of  the  cold  front) ,  sometimes 
appearing  as  an  extension  from  the  duct  region  associated  with 
the  gulf  stream.  During  the  Southern  Hemisphere  winter,  ducting 
is  frequently  identified  in  the  South  Pacific  in  the  area  south 
west  of  the  Chilean  coast  near  Antarctica.  This  is  an  area  with 


a  high  frequency  of  cyclogenesis.  These  storms  are  associated 
with  the  advection  of  extremely  cold  air  off  of  the  Antarctic 
continent  over  relatively  warm  ocean  surfaces.  Thus,  these 
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storms  may  activate  the  model  cold  air  advection  parameterization 
in  which  the  diagnosed  air/sea  temperature  difference  is  replaced 
by  a  prescribed  cooling  rate  based  on  the  geostrophic  advection 
of  air/sea  temperature  difference.  Large  near  surface  gradients 
in  temperature  contribute  to  near  surface  N  gradients 
characteristic  of  surface  evaporative  ducts. 

Examples  of  the  general  features  discussed  above  appear  in 
Figures  2-13..  Figures  2-9  show  a  series  of  \2  hour  sea  level 
pressure  (SLP)  forecasts  with  corresponding  12  hour  forecasts 
of  SED  beginning  at  00Z  on  29  January  and  ending  at  12Z  on  31 
January  1987.  An  intense  storm  is  located  in  the  Gulf  of  Alaska 
throughout  the  period.  There  is  no  indication  of  any  ducting 
anywhere  near  the  Gulf.  In  Figure  2,  note  the  storm  in  the  North 
Atlantic  near  40  degrees  N  and  40  degrees  W.  In  Figure  3,  note 
the  extension  of  the  duct  region  near  the  east  coast  of  North 
America  (near  the  Gulf  Stream)  towards  the  location  of  the  storm. 
An  example  of  a  ducting  area  on  the  southeast  side  of  a  low  in 
the  North  Atlantic  is  shown  in  Figures  10  and  11  for  00Z  on  21 
January  1987.  In  this  situation,  the  duct  area  appears  to  be 
nearly  coincident  with  the  area  of  maximum  pressure  gradient. 
Figures  14-16  show  a  series  of  forecasts  of  SED  all  valid  at 
00Z  on  26  June  1986.  The  analysis  is  shown  in  Figure  12  and 
the  corresponding  SLP  analysis  in  Figure  13.  This  series  has 
several  interesting  features,  including  responses  to  storms  near 
Antarctica.  The  size  and  depth  of  the  duct  area  between  the 
southwest  coast  of  Peru  and  Antarctica  are  quite  consistent 


Figure  10.  12  hour  forecast  of  the  sea  level  pressure  field  valid  at 

00Z  on  21  January  1987. 


Figure  11.  12  hour  forecast  of  the  radar  evaporative  duct  height  field 

valid  at  00Z  on  21  January  1987. 


L_ 


^ 

v;  /?/.v 


V_  AV. 


23 


Figure  12.  Analysis  of  the  radar  evaporative  duct  height  field  valid 
at  00Z  on  26  June  1986. 


June  1986. 
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ure  15.  24  hour  forecast  of  the  radar  evaporative  duct  height  field 

alid  at  00Z  on  26  June  1986. 


throughout  the  period  and  agree  well  with  the  analysis,  although 
the  shape  and  exact  location  are  more  variable.  Figure  13  shows 
an  intense  storm  due  south  of  Peru.  The  duct  area  off  the  coast 
of  Peru  appears  to  be  coincident  with  the  area  of  maximum  winds. 
Note  also  the  storms  south  of  Australia  and  in  the  south  Atlantic 
ocean  in  Figure  13.  Figures  12  and  14-16  show  that,  while  there 
is  ducting  south  of  Australia,  the  south  Atlantic  is  duct-free 
throughout  the  period. 

3.3  Tropics  and  Subtropics 

The  response  of  the  SED  product  in  the  tropics  and  sub¬ 
tropics  is  typified  by  few  forecasts  of  ducting  in  the  central 
Indian  ocean,  central  equatorial  Pacific,  and  west  equatorial 
Atlantic  with  somewhat  better  performance  in  the  eastern  Indian 
ocean  and  near  the  Asian  continent.  There  is  no  significant 
diurnal  difference. 

In  general,  the  ducting  situation  in  the  model  is  quite 
persistent  with  certain  geographical  areas  favored  for  SED 
formation  and  certain  other  areas  in  which  ducts  rarely  form. 

With  the  exception  of  the  equatorial  Pacific,  all  of  the  favored 
areas  are  near  continents  -  east  of  Australia  and  Africa  and 
in  the  South  China  Sea  and  Bay  of  Bengal .  In  the  equatorial 
Pacific,  there  is  a  bias  toward  more  extensive  duct-filled  areas 
for  longer  forecast  periods.  The  persistence  of  the  ducting 
situation  in  the  tropics  is  not  surprising;  the  weather  in  the 
tropics  is  also  quite  persistent.  It  is  important  to  note  that 


there  is  almost  no  ducting  activity  forecast  or  analyzed  in  any 
of  these  summer  cases  off  the  coast  of  California,  which  is  also 
an  area  favored  for  SED  formation. 


Some  examples  of  the  features  discussed  above  from  the 
summer  of  1986  appear  in  Figures  12-21.  In  comparing  Figures 
14-16  (successive  forecasts)  with  Figure  12  (the  analysis  from 
00Z  26  June  1986),  several  important  features  are  evident.  For 
example,  the  extensive  region  of  duct  heights  greater  than  16  m 
extending  across  the  Indian  ocean  and  surrounding  Australia  in 
the  analysis  is  present  in  the  12  hour  forecast,  present  only 
east  of  Australia  in  the  24  hour  forecast,  and  completely  absent 
in  the  36  hour  forecast.  Although  forecast  duct  heights  greater 
than  16  m  are  present  in  the  central  equatorial  Pacific  through¬ 
out  the  period,  the  exact  location  is  quite  variable .  Some 
features  which  are  better  forecast  are  the  area  off  the  east 
coast  of  Africa  from  about  5  degrees  S  to  about  25  degrees  S 
and  the  area  in  the  South  China  Sea  and  Bay  of  Bengal .  The  area 
of  ducting  in  the  equatorial  Atlantic  off  the  coast  of  Brazil 
appears  consistently  in  all  forecasts  but  is  not  realized;  i.e., 
it  is  not  in  the  analysis . 

Many  of  the  same  features  are  present  in  a  series  of 
forecast  valid  at  12Z  on  24  June  1986  (not  shown)  and  in  a 
similar  series  for  12Z  on  30  June  1986  (Figures  17-21)  .  In  the 
30  June  series,  ducting  across  the  Indian  ocean  is  not  present 
although  the  area  of  ducting  off  the  east  coast  of  Africa  is. 
Ducting  east  of  Australia  is  also  present.  Ducting  in  the 
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central  equatorial  Pacific  is  somewhat  more  extensive  them  in 
the  26  June  series  in  both  the  analysis  and  the  forecasts  through 
24  hours  (at  hour  36,  the  area  is  not  as  extensive  as  in  the  26 
June  series) .  Also,  the  area  extends  further  to  the  north  (north 
of  15  degrees  N) .  The  ducting  area  in  the  South  China  Sea  is  not 
as  apparent  in  the  analysis,  however,  the  area  in  the  Bay  of 
Bengal  is  present  in  the  30  June  series.  Even  the  nonrealizable 
feature  in  the  equatorial  Atlantic  is  present  in  this  series. 

4.  CONCLUSION  AND  RECOMMENDATIONS 

Surface  evaporative  ducts  have  considerable  tactical  signif¬ 
icance  in  terms  of  the  threat  of  detection  by  hostile  units  due 
to  extended  range,  loss  of  the  ability  to  detect  hostile  aircraft 
due  to  "holes",  and  danger  of  enhanced  clutter.  The  only 
routinely  available  forecast  aid  for  SED  prediction  is  the  SED 
product  from  the  global  surface  contact  layer  interface  model . 

The  GSCLI  model  is  a  generalized  similarity  theory  model  for  the 
entire  PBL  in  which  the  NOGAPS  low  level  pressure  gradient,  sea 
surface  temperature,  and  850  mb  temperature  are  used.  The  method 
of  stability  determination,  constraints  placed  on  important 
parameters,  and  "engineering  fixes"  employed  in  the  GSCLI 
warrant  an  examination  of  the  quality  of  the  SED  product. 

This  e'  aluation  of  the  SED  product  indicates  that  there 
is  a  considerable  geographical  bias  in  the  quality  of  the  SED 
forecasts.  In  the  midlatitudes,  the  North  Atlantic  and  South 
Pacific  regions  are  well  supported  by  the  SED  product  while,  in 
the  North  Pacific  and  South  Atlantic,  evaporative  ducting  is 
rarely  forecast,  regardless  of  the  synoptic  situation.  In  the 
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tropics  and  subtropics,  persistent  ducting  conditions  near 
continents  are  well  forecast  (except  near  Brazil,  where  forecast 
ducting  almost  never  materializes)  while  ducting  over  open 
equatorial  waters  is  poorly  forecast . 

It  should  be  noted  that,  since  the  GSCLI  is  a  purely 
diagnostic  system,  some  of  the  apparent  deficiencies  in  forecast 
fields  may  be  reflections  of  deficiencies  in  the  NOGAPS  forecast 
fields  used  by  the  GSCLI .  There  is  very  little  data  in  the 
tropics  and  most  of  it  covers  tropical  continental  areas.  In 
the  midlatitudes,  extensive  analyses  of  NOGAPS  forecast  fields 
performed  at  the  FNOC  have  shown  a  consistent  cold  bias  in  near¬ 
surface  air  temperatures  in  the  North  Pacific  ocean  at  all 
forecast  times  to  72  hours  (Quality  Control  Division,  1986) . 

This  may  play  a  role  in  the  geographical  bias  of  the  A90  field 
in  the  North  Pacific.  No  analyses  have  been  performed  on  the 
Southern  Hemisphere  fields,  but  it  is  conceivable  that  similar 
behavior  may  be  responsible  for  A90  performance  in  the  South 
Atlantic . 

It  should  also  be  noted  that,  since  these  analyses  of  the 
surface  duct  height  product  were  completed,  a  new  global  analysis 
and  forecast  system  (NOGAPS  3.0)  has  become  operational  at  the 
FLENUMOCEANCEN .  Therefore,  some  of  the  conclusions  discussed 
herein  regarding  geographic  variability  may  no  longer  be  valid. 

An  analysis  similar  to  this  one  could  be  performed  in  order  to 
assess  the  impact  of  NOGAPS  3.0  on  the  surface  evaporative  duct 
height  product . 


Based  on  the  arguments  presented  in  Appendices  A  and  B,  it 
is  apparent  that  the  expression  used  in  the  GSCLI  model  for  the 
surface  evaporative  duct  height  possesses  a  serious  shortcoming. 
When  the  ambient  conditions  are  such  that  dN/dz  never  approaches 
the  critical  value,  the  GSCLI  expression  falsely  produces 
reasonable  appearing  duct  heights.  This  seems  to  be  undesirable. 
It  is  recommended,  therefore,  that  the  similarity  theory  approach 
for  potential  refractivity  presented  in  Appendix  B  be  adopted 
under  stable  conditions  for  the  operational  surface  evaporative 
duct  height  product . 
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APPENDIX  A 


DERIVATION  OF  SURFACE  EVAPORATIVE  DUCT  HEIGHT 


The  expression  for  6  In  the  GSCLI  model  can  be  obtained  by 


taking  the  vertical  derivative  of  the  equation  defining  the  index 


of  refraction  N: 


N  =  APT  1  +  APq ( B/ 1 ) T~  2 


dN  3N  dT  3N  dq  3N  dP 

—  =  —  —  +  —  —  + - 

dz  3T  dz  3q  dz  3P  dz 


Evaluating  the  partial  derivatives  in  (2)  gives 


=  -APT  2  -  2APq(B/€)T  3 


- - 2  C 1  +  2q(B/eT)] 

T 


3N  AP 


-  =  -=■  (B/e) 

3q  x 


3N  A  AqB 


3P  T  sT 


=  -  [1  +  q(B/eT) ] 
T 


thus , 


dP/dz  =  -p g  =  -Pg/RT 


3N  dP  AP 

- - - ~  [1  +  q(B/-:T)  ]  (g/R) 


Using  (3),  (4),  and  (5),  we  can  re-write  (2)  as 


dN 

AP 

dT 

dq 

—  = 

[ (1  +  2q ( B/eT ) )  — 

-  (B/e)  — 

dz 

T2 

dz 

dz 

+  (1  +  q(B/sT) ) g/R ]  (6) 

The  remaining  substantial  derivatives  in  (6)  can  be  evaluated 
from  surface  boundary  layer  similarity  expressions: 


dT 

T* 

=  —  +(z/L) 

(7) 

dz 

KZ 

dq 

<3* 

=  —  *(z/L) 

(8) 

dz 

KZ 

Substituting  (7)  and  (8)  into  (6)  gives 


dN  AP  *(z/L) 

-  - - r-  [{(1  +  2q(B/eT))T*  -  (B/e)q*}  - 

dz  kz 

+  (1  +  q ( B/eT ) ) g/R ]  (9) 


At  the  duct  height  the  gradient  in  N  is  equal  to  the  critical 
value  b2  : 


AP  * ( o /L ) 

b,  - - -  [{(1  +  2q(B/'=T)  )T*  -  (B/e)q*>  - 

+  (1  +  q ( B/^T ) ) g/R ]  (10) 


At  this  point,  two  simplifications  are  made  in  the  derivation. 
First,  in  the  stability  function  *(<f/L),  6  is  replaced  by  Z,  a 
reference  level.  Retaining  S  in  4>(^/L)  under  unstable 
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conditions  would  require  an  iterative  solution.  Second,  terms  of 
order  qB/*T  were  neglected  -  this  is  difficult  to  justify;  for  a 
specific  humidity  of  10  g/Kg  and  an  air  temperature  of  300  K, 
qB/eT  =  0.26,  2qB/*T  =  0.52  which  is  certainly  not  small  in 
comparison  to  1  and,  in  the  tropics,  2qB/sT  may  be  ~0(1).  We 
will  proceed  with  the  derivation  of  <S  and  then  derive  an  alternate 
expression  in  which  neither  of  the  above  simplifications  are 
made.  With  the  above  simplifications,  equation  (10)  becomes 


AP 


* (Z/L) 


b„  - - r  HT*  -  (B/*)q*> 


K6 


+  g/R] 


(11) 


K<fb2  —  +  g/R  K $  =  (B/e)q*  <MZ/L)  -  T*  <t>(Z/L) 
AP 


(B/e)q*  (  Z/L  )  -  T*  *<  Z/L  ) 


0  = 


k( (b£T  /AP)  +  g/R) 


(12) 


Equation  (12)  is  the  expression  used  in  the  GSCLI  model. 
Returning  now  to  the  full  expression  (10),  we  can  derive  an 
alternate  expression  for  <f : 


<f  = 


-  [{(1  +  2q{  B/'iT  )  )  T*  -  (  B/e  )  q*  }  <M<f/L)] 


K[(b„T  /AP)  -  (1  +  q ( B/eT ) ) g/R ) 


(13) 


Under  stable  conditions, 


<Mz/L) 


0.74  +  4. 7 (z/L) 


(14) 
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Substituting  (14)  into  (13)  gives 

-  [{(1  +  2q(  B/sT  )  )T*  -  (B/s  )q*}  (  0 . 74  +  4.7U/L))] 

4  - - 

D 

where  D  is  the  denominator  of  (13);  i.e., 

D  =  K[(b2T2/AP)  -  (1  +  q ( B/sT ) ) g/R ] 

4  =  -  1 0 . 74 { ( 1  +  2q(B/eT))T*  -  (B/€)q*> 

+  4 . 7  ( <f/L )  {  ( 1  +  2q(B/eT))T*  -  (B/e)q*J)  D_1 

-0 . 74  { ( 1  +  2q  ( B/'sT )  )  T*  -  ( B/ >= )  q*  } 

4  =  -  (15) 

D  +  4.7/LU  +  2q(B/eT))T*  -  (B/e)q* 

Under  unstable  conditions, 

*(z/L)  =  0.74(1  -  9(z/L))~1/2  (16) 

Substituting  (16)  into  (13)  gives 

4  =  -  [((I  +  2q(B/eT))T*  -  (B/e)q*} 
x  (0.74(1  -  9U/L)  )"1/2]  D_1 

<f  ( 1  -  9U/L))~1/2  =  -  0 . 74  [  ( 1  +  2q  (  B/sT  )  )  T* 

-  ( B/e ) q* ]  D_1  (17) 


Equations  (15)  and  (17)  are  considerably  more  complicated  than 
(12).  In  order  to  assess  the  significance  of  the  differences 
between  (12)  and  (15)  and  (17),  a  program  was  written  which  eval¬ 
uates  (12)  and  either  (15)  or  (17),  depending  on  atmospheric 


stability.  Results  produced  by  the  program  revealed  that,  in 
near-neutral  conditions,  the  two  expressions  produced  very  nearly 
identical  results.  For  increasing  stability,  the  full  expression 
(15)  gave  duct  heights  substantially  larger  than  the  expression 
(12)  used  in  the  GSCLI  model.  Under  unstable  conditions,  the 
Newton-Raphson  technique  was  used  to  solve  (17).  Solutions  were 
obtained  in  3-5  iterations  with  an  initial  estimate  of  10  m. 
Equations  (17)  and  (12)  produced  similar  results.  Results  also 
revealed  an  important  characteristic  of  equation  (12).  An  as¬ 
sumption  inherent  in  (10)  is  that,  somewhere  in  the  surface  lay¬ 
er,  dN/dz  =  b2;i.e.,  that  a  surface  evaporative  duct  exists. 

When  this  is  not  the  case,  (12)  still  produces  values  that  appear 
to  be  resonable  for  duct  height.  In  reality,  when  dN/dz  <<  [b£l, 
no  duct  exists  and  when  dN/dz  >>  [b£],  the  surface  evaporative 
duct  height  is  greater  than  the  depth  of  the  surface  layer  and 
there  is  strong  trapping  extending  up  the  elevation  at  which 
dN/dz  <  [b£l.  Equations  (15)  and  (17)  produce  very  shallow  duct 
heights  when  dN/dz  >>  [b£]  and  large  negative  duct  heights  when 
dN/dz  <<  [b„].  The  fact  that  equation  (12)  produces  false  indic¬ 
ations  of  reasonable  duct  thickness  in  these  situations  seems 
undesirable . 


APPENDIX  B 


APPLICATION  OF  SURFACE  SIMILARITY  THEORY  TO 
POTENTIAL  REFRACT IVITY 

In  this  Appendix  we  first  demonstrate  that  potential 
refractivity  cam  be  scaled  by  the  same  type  of  surface  layer 
similarity  expressions  as  commonly  used  for  temperature  and 
moisture  (ordinary  refractivity  cannot  be  treated  this  way, 
however) .  Then,  using  this  similarity  expression  for  potential 
refractivity,  we  derive  expressions  for  evaporative  duct  height 
as  a  function  of  stability.  This  method  of  computing  duct 
height  is  then  commpared  with  others  discussed  in  this  report. 

Similarity  theory  as  applied  to  mean  and  turbulence 
variaibles  shows  that  many  parameters  follow  simple  scaling  laws 
in  the  atmosheric  surface  layer.  The  functional  forms  of  these 
scaling  laws  have  been  fairly  well  established  by  observational 
studies.  For  example,  in  the  surface  layer,  nondimensional 
temperature  and  moisture  gradients  follow  the  scaling  laws: 

K  z  9T 

—  —  =  (Z/L)  ,  (1) 

T*  3z 

Kz  9q 

-  —  “  <P  Q  <Z/L>  '  (2) 

q*  3z 

where  T*  and  q*  are  scaling  parameters  and  L  is  the  Monin- 
Obukhov  length.  Generally  the  universal  functions  for  heat  and 
moisture  are  taken  to  be  equal,  i.e.,  ^  =  4*^  =  4>  . 
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It  is  often  assumed  (without  proof)  that  refractivity  cam 
be  scaled  in  the  sarnie  manner,  so  that 


k  z  8  N 

-  -  =  4,  (Z/L)  , 

N*  3z 


but  we  will  see  that  this  is  wrong.  The  expression  for 
refractivity  in  the  microwave  region  of  the  spectrum  is 


A  r 

N  =  -  IP  + 


A  quantity  called  potential  refractivity,  X,  can  be  defined  such 
that  it  does  not  vary  with  pressure : 


A  r  BP0q, 

"  -  Po  +  - 

T  l  £T  J 


where  PQ  =  1000  mb.  A  total  differential  of  X  is 

dX  =  —  dT  +  —  dq  . 

3  T  3  q 


Replacing  the  differentials  in  (6)  with  fluctuation  notation 
(a  common  trick  in  turbulence  derivations)  yields: 


•  9X  , 

X  =  - —  T'  +  - —  q  . 


If  we  multiply  (7)  by  -w'  divide  by  U*,  and  average: 


w' X'  3X  w'T'  3X  w'q' 

-  =  (_)  (-  - )  +  (_)  (-  - ) 

U*  3T  U*  3q  U* 


or,  using  the  standard  notation  of  scaling  parameters, 

ax  m  ax 

X*  =  -  T*  +  -  q*  . 

a  t  a  q 


Returning  now  to  (6)  and  writing  x  as  a  total  vertical 
derivative. 


dX  3X  dT  3X  dq 

—  =  -  —  + - 

dz  .  3T  dz  3q  dz 


and  substituting  (1)  and  (2)  for  the  temperature  and  moisture 
gradients  gives : 


dX  f  ax  a X  1  <P  (z/L) 

—  =  (— )T*  +  (— )qJ  - 

dz  l  3t  3q  J  kz 


Now,  using  (9)  we  can  rewrite  (11)  as: 


dx  <Mz/L) 

—  -  X*  - 

dz  k  z 


<z  dX 

- - =  <M  z/L) 

X*  dz 


Thus,  we  have  proved  that  the  potential  refractivity  obeys 
surface  layer  similarity. 

Note,  however,  that  had  we  attempted  this  derivation  with 
N  rather  than  Xr  we  would  have  found  in  place  of  (8) : 


w' N'  3N  w'T'  3N  w' q'  3N  w'p' 

-  =  ( - )  (-  — )  +  (_)  (-  - )  +  (_)  (-  - )  ( 

U*  aT  u*  3q  u*  3P  u* 


and  we  could  not  have  easily  treated  the  last  term  in  this 
expression.  Thus,  it  is  proper  to  treat  x(T*q)  a3  a  similarity 
variable,  but  not  N(T,q,p). 

We  know  that  the  critical  refractivity  gradient  for 
trapping, 

dN 

dz  c, 

is  -0.157  m-^.  Now  we  need  to  know  what  this  corresponds  to 
in  terms  of 


First,  by  examination  of  (4)  and  (5)  we  note  that 


3N 

9X 

3N 

9x 

-  = 

=  - 

and  —  *  — 

3T 

3T 

9q 

aq 

Taking  a  total  derivative  of  (4)  gives 


so  that. 


dN 

3N 

dT 

3N 

dq 

3N 

dp 

_  — 

-  +  - 

-  +  - - 

dz 

9T 

dz 

9q 

dz 

9P 

dz 

dN 

3X 

dT 

9x 

dq 

3N 

dp 

_  —  = 

-  +  - 

-  +  - - 

dz 

3T 

dz 

9q 

dz 

9P 

dz 

^^^vjvvatrv.v.v'^  'rw.v.w  v v v  v**", 


dp 


Using  the  hydrostatic  equation,  — 

dz 


pg 


RdT 


gives 


dx  dN  pg  3N 


(18) 


dz  c  dz 
From  (4)  we  find 


c  —  c  RdT  3P 


3N  A  Bq 

—  -  -  [1  +  — ]  , 

3p  T  £T 


(19) 


so  that. 


dN 


( — )  =  (_)  + 
dz  „ 


pgA 


Bq 


[1  +  — 1  • 


(20) 


Thus,  the  critical  x  gradient  is  not  a  constant,  but  rather  a 
function  of  p,  T,  and  q.  We  can  estimate  a  typical  magnitude  of 
this  last  term  by  selecting  p  =  1000  mb,  T  =  300°k,  and  q  =  10-2. 
Then 


pg- 


RdT‘ 


Bg- 

1+  —  |~  0.037. 
T. 


In 

-  ~  0 

rt 


Thus, 


d  X 


( — )  =  -0.157  +  0.037  =  -0.120  m 

dz  „ 


-1 


The  IREPS  model  uses  a  constant  value  of  -0.125  m-^  for  the 
critical  gradient  of  refractivity . 


We  purposely  have  not  written  the  hydrostatic  equation  in  terms 
of  density  because  the  coefficients  A  and  B  are  dimensional  and 
the  values  we  are  using  are  designed  to  be  used  in  conjunction 
with  pressure  in  ■  . 


48 


39 


& 


:fll 


Q 


Wl 


A 


> 


.y 


*v_ 


'  V 


1 
i 


5 
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if.  A  A  j*./ .  t  .V.  <.  f. 


dX 

If  we  let  Y  =  ( — )  ,  and  treat  this  as  a  constant,  we  may 

„  dz  c 

write  (12)  as: 


<5  r 

—  (Y)  =  <M-) 

X*  L 

where  6  is  the  duct  height  at  which  X  reaches  its  critical 
gradient.  For  4>  (z/L)  we  take  [Businger,  1973] 


(21) 


4>(Z/L)  = 


.74  +  4.7 (z/L) 


L  >  (neutral,  stable) 


.74(1  -  9  z/L)-1/2  L  <  0  (unstable 


(22) 


Substituting  (22)  into  (21)  we  find,  for  stable  conditions: 

7x*r1 

-  (23) 

Yl  / 


/KY  4.7^ 

|  .  '74X-  | 

f!  -i 

L  ' 

<  Y 

^  < 

[At  this  point  we  only  note  in  passing  that  the  expression 

for  <5  can  go  negative  under  certain  conditions.] 

For  unstable  conditions  we  find, 

6  1/2  • 74  X  * 

6(1  -  9  -  )  =  - 

L  <  Y 


(24) 


For  ease  of  use  with  a  Newton-Raphson  technique,  we  square  this 
equation  and  rearrange  as : 


F  (6 )  =  6 2  ( 1  -  9  -)  -  ( 


.  74X, 


->  =  0 


(25) 


L  <  Y 

The  derivative  of  this  equation  is 


F' (6)  =  26 


-3f-) 

\T.  / 


(26) 


49 


and  thus  we  can  solve  for  <5  by  the  Newton-Raphson  iterative 
method : 


» 

*  * 
►  s 

’.s' 

£ 


5 (m+1) 


6(m) 


F(fi(m)) 
F'  (6(m) ) 


(27) 


where  is  the  (m+1)  approximation  to  6.  This  converges 

rapidly.  Thus,  we  can  solve  for  duct  height  by  use  of  either 
(23)  or  (24)  provided  we  have  X*.  And  X*  can  be  found  by  first 
noting  (analogous  to  T  and  q)  we  can  write  the  integrated  form 
of  (12)  as 

•  74X* 

X  (z)  -  X(zQ)  - -  (In  z/zQ  -  ip  (z/L)  ]  (28) 

< 


u 


or,  rewriting  as, 


K (X(Z)  -  X(zo) ] 

X*  =  - 

.74  [In  z/zQ  -  ip  (z/L) 


(29) 


Here  ip(z/L)  are  integrated  forms  of  0(z/L)  and  have  been  found 
to  be 


ip(  z/L) 


4.7 
.74 
2  In 


z/L 

(l+y) 

2 


,  L  >  0 

,  L  <  0 


(30) 


where,  y  =  (1-9Z/L)  ^  . 

Now  we  summarize  the  steps  for  computing  the  evaporation 
duct  height  <5  based  on  a  similarity  expression  for  potential 
refractivity,  X.  we  assume  that  we  already  know  values  of  T 
and  q  at  some  reference  height,  <r,  and  we  know  the  sea  surface 


a 


v  -*i 


! 


>3 

I 


:-3> 

o 

k 

a? 


temperature.  Further,  we  assume  that  we  have  already  computed 
L  and  zQ  from  expressions  that  we  do  not  discuss  here. 

1.  Compute  X(zr)  and  X(zQ)  using  (5). 

2.  Compute  X*  from  (29)  with  z-*-zr. 

3.  Depending  on  the  sign  of  L,  compute  6  either  from  (23) 
or  (24)  . 

Of  the  different  methods  discussed  in  this  report  for 
finding  5,  the  method  outlined  above  appears  to  be  the  most 
accurate  and  to  require  relatively  little  computing  power. 
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